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ABSTRACT 

We demonstrate a successful strategy for identifying extremely metal poor galaxies. Our preliminary 
survey of 24 candidates contains 10 metal poor galaxies of which 4 have log(0/II) + 12< 7.65, some 
of the lowest metallicity blue compact galaxies known to date. Interestingly, our sample of metal 
poor galaxies have systematically lower metallicity for their luminosity than comparable samples of 
blue compact galaxies, dirrs, and normal star-forming galaxies. Our metal poor galaxies share very 
similar properties, however, with the host galaxies of nearby long-duration gamma-ray bursts (GRBs), 
including similar metallicity, stellar ages, and star formation rates. We use II/3 to measure the number 
of OB stars present in our galaxies and estimate a core-collapse supernova rate of ^^10^^ yr^^. A 
larger sample of metal poor galaxies may provide new clues into the environment where GRBs form 
and may provide a list of potential GRB hosts. 

Subject headings: galaxies: abundances — galaxies: starburst — gamma rays: bursts 



> 

o 
o 

o 
o 



X 



1. INTRODUCTION 

Metal poor galaxies are the key to understanding star 
formation and gas enrichment in a nearly pristine in- 
terstellar medium, and may provide a template for un- 
derstanding the formation of the first stars. Extremely 
metal poor galaxies (XMPGs) are extremely rare: fewer 
than 1% of dwarf galaxies are XMPGs, with a gas-phase 
oxygen abundance log(0/H) + 12< 7.65 (Kunth & Ostlin 
2000; Kniazev et al. 2003). Known XMPGs are mostly 
gas-rich, blue compact galaxies with spectra dominated 
by emission lines. The first surveys to search for XMPGs 
were objective prism surveys (MacAlpine et al. 1977; 
Kunth et al. 1981; Terlevich et al. 1991). Abundance 
studies of these surveys revealed up to a dozen XMPGs 
(Kunth & Sargent 1983; Campbell et al. 1986; Masegosa 
et al. 1994), but none so extreme as I Zw 18 (Searle 
& Sargent 1972). I Zw 18, now along with SBS 0335- 
052W and DDO 68, are the most metal poor local galax- 
ies known, with log(0/H) -|- 12 ranging 7.12 - 7.17 (Izo- 
tov et al. 2005; Izotov & Thuan 2007). Papaderos et al. 
(2006) found 2 new XMPGs in the 2dF survey. Among 
the 1,000,000 spectra released by the Sloan Digital Sky 
Survey (SDSS), there are 19 identified XMPGs (Kniazev 
et al. 2003; Izotov et al. 2004, 2006a; Izotov & Thuan 
2007), emphasizing the rarity of such objects. 

Kewley et al. (2007) recently announced the 
serendipitous discovery of a new XMPG, SDSS 
J080840.85-fl72856.48 (hereafter SDSS J0809-I-1729). 
The object is a stellar point source in the SDSS cata- 
log. It was observed as part of the Brown et al. (2006a) 
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hypervelocity star survey on the basis of its stellar B- 
type photometric colors. Spectroscopy reveals that the 
object is not a star but rather a compact blue galaxy 
at cz = 13232 km s^^. Our re-analysis shows that the 
galaxy has log(0/II) + 12= 7.48 ± 0.1. Interestingly, the 
observed electron density, star formation rate, and to- 
tal luminosity of this XMPG are remarkably similar to 
nearby GRB host galaxies (Kewley et al. 2007). 

Nearby z < 0.2 long duration GRBs are observed 
in metal poor galaxies (Prochaska et al. 2004; SoUer- 
man et al. 2005; Fruchter et al. 2006; Stanek et al. 
2006; Wolf & Podsiadlowski 2007; Wiersema et al. 2007; 
Margutti et al. 2007). Fruchter et al. (2006) argue that 
the link between metal poor galaxies and GRBs origi- 
nates in the atmospheres of massive stars. Massive metal 
poor stars lack the opacity to support significant stellar 
winds, and thus can produce the anomalous Type Ic su- 
pernovae associated with nearby GRBs. Berger et al. 
(2007) argue that GRBs are linked to young starburst 
populations, which at low redshift happen to be found 
predominantly in low mass, metal-poor galaxies. At 
large redshift z > 0.2, the metallicity of GRB hosts is 
much harder to determine (Prochaska 2006). At least 
one GRB host galaxy is an extremely red and proba- 
bly metal rich object (Berger et al. 2007), and other 
GRB hosts 0.2 < z < 1 appear consistent with normal 
metallicity-luminosity relations (Wolf & Podsiadlowski 
2007; Margutti et al. 2007). GRBs clearly occur in dif- 
ferent types of host galaxies; they are not tied exclusively 
to the most metal poor galaxies. Yet finding and study- 
ing metal poor galaxies with properties similar to nearby 
GRB host galaxies may yield critical clues about the en- 
vironment where nearby GRBs form. 



2 



Brown, Kewley & Geller 



2.5x10^'!' 




3000 



4000 



5000 6000 
X [A] 



7000 



8000 



Fig. 1.— Observed spectrum of the XMPG SDSS J0809+1729 
(Kewley et al. 2007) plotted against the total system throughput 
of the SDSS u' g' r' i' filter bandpasses (dotted lines). The strong 
emission line contribution to the galaxy's broadband magnitudes 
results in unusual colors that vary with redshift. 



Inspired by the discovery of SDSS J0809+1729, we de- 
signed a survey to find metal poor galaxies. We use 
a technique similar to photometric rcdshifts to identify 
metal poor galaxies in the SDSS galaxy catalog. Wc test 
this technique and target g' ^ 20 galaxies with very blue 
colors, a region of parameter space not well probed by 
previous surveys. Our strategy uncovers 10 metal poor 
galaxies from a sample of 24 candidates, 4 of which are 
new XMPGs. 

In §2 we present our technique to find new metal poor 
galaxies and discuss the efficacy of our survey. In §3 we 
describe the properties of the entire set of metal poor 
galaxies, and compare the galaxies with samples of blue 
compact galaxies and nearby GRB hosts. In §4 we es- 
timate the expected core-collapse supernova rate in our 
galaxies. We conclude in §5. 

2. DATA 

2.1. Technique to Find Extremely Metal Poor Galaxies 

Known XMPGs outside the Local Group are starburst 

galaxies characterized by very low internal extinction, 
high ionization parameter, and large gas-phase electron 
density (e.g. Kniazev ct al. 2003; Izotov et al. 2005, 
2006a; Papaderos et al. 2006; Kewley et al. 2007; Izo- 
tov & Thuan 2007). As a result, XMPGs have steep 
blue continua, large [Olll]/[Oll] ratios, and strong hy- 
drogen Balmer emission lines. These characteristics sig- 
nificantly affect the broadband colors of XMPGs. For 
example. Figure 1 shows the observed spectrum of SDSS 
J0809+1729 plotted against the total system throughput 
of the SDSS u'g'r'i' filter bandpasses^. The contribution 
of emission lines to the broadband magnitudes of SDSS 
J0809-hl729 is approximately 3%, 32%, 6%, and 12% at 
u', g', r' , and i' respectively. 

The strong emission line contribution to XMPGs' 
broadband magnitudes results in unusual colors that 
change with redshift. For example, at lower redshift, 
SDSS J0809+1729 docs not satisfy the B-star criteria in 
Brown et al. (2006b) because Ha drops into the r' band 
and produces a much redder {g' — r') color. Conversely, 
at higher redshift, H/3 and O III move from 9' into r' and 
also produce a much redder {g' — r') color (see Figure 
1). We quantify these effects using the RVSAO pack- 

^ http://www.sdss.org/dr5/instruments/imager/index.html 
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Fig. 2. — Target selection for our MMT survey for XMPGs. 
a) Color-color plot showing all SDSS DR4 galaxies (dots) and the 
color-redshift track for the XMPG J0809-I-1729 (solid line), b) - d) 
Color-redshift plots showing SDSS objects (dots) that simultane- 
ously match all three color-redshift tracks (solid lines) within 0.1 
mag. After removing saturated stars and nearby Hli regions by vi- 
sual inspection, we are left with a sample of 38 XMPG candidates 
(x's). Four new XMPGs discovered by this survey are indicated by 
stars. 



age (Kurtz & Mink 1998) to shift the XMPG spectra to 
different redshifts. We predict the galaxy's colors at dif- 
ferent redshifts by adding the relative change in color to 
the observed photometry. In effect, we are calculating 
^-corrections for the XMPG. 

Figure 2 plots the resulting color-redshift tracks for 
SDSS J0809-hl729 (sohd lines). Over the range 0.0 < 
z < 0.10, the XMPG's («/ - ,g')o and {g' - r')o colors 
vary by ^ 0.3 mag and the (r' — i')o color varies by 
^1 mag. The subscript indicates colors corrected for 
Galactic extinction following Schlegel et al. (1998). We 
note that the Johnson/Cousins passbands are much less 
sensitive than their SDSS equivalents because the strong 
H/3/[0 111] and Ha emission lines remain in their respec- 
tive V and R passbands for 0.0 < z < 0.10. 

Knowing how the broadband colors of an XMPG 
change with redshift, we can search for new XMPGs 
at other redshifts in the SDSS photometric catalog. 
Unfortunately, searching for new XMPGs like SDSS 
J0809-f 1729 in the stellar catalog is not feasible because 
of immense contamination from white dwarfs with simi- 
lar colors. Instead, we search the SDSS galaxy catalog. 

We use the color-redshift track of our newly discovered 
XMPG, SDSS J0809-I-1729, as the basis for the sample of 
XMPGs published here. We begin by selecting all galax- 
ies in SDSS DR4 (Adelman-McCarthy et al. 2006) with 
g' < 20.5. Figure 2a plots these galaxies with similar 
{u' - g')o and {g' - r')o to SDSS J0809-I-1729. Only by 
combining two or more colors can we meaningfully select 
XMPG candidates. We find ~ 10^ galaxies simultane- 
ously within 0.1 mag of any pair of color-redshift tracks, 
and a mere 107 galaxies (small squares. Figure 2b-d) si- 
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multaneously within 0.1 mag of all three color-redshift 
tracks. We visually inspect the objects and find that 
many are near saturated stars or that they are HII re- 
gions in nearby galaxies. After eliminating the unwanted 
objects, we are left with a sample of 38 photometrically- 
selected XMPG candidates (marked by x's, Figure 2b-d). 
Our survey is based on this sample of 38 XMPG candi- 
dates. 

We use a restrictive color selection as the first demon- 
stration of our technique. If our XMPG selection strat- 
egy is successful, we can easily broaden the search pa- 
rameters to identify many more faint XMPG candidates. 

2.2. Observations 

We obtained spectroscopy of the 24 candidates avail- 
able on the nights of 2006 May 25-27 and 2006 June 19- 
20. Table 1 lists the 24 candidates. Observations were 
obtained with the 6.5m MMT telescope and the Blue 
Channel spectrograph. We operated the spectrograph 
with the 300 line mm~^ grating and a 1" slit. These set- 
tings provide a wavelength coverage of 3400 A to 8600 A 
and a spectral resolution of 6.2 A. Exposure times were 
typically 30 minutes. We obtained comparison lamp ex- 
posures after every exposure. 

We reduce the spectra using standard IRAF^ spec- 
tral reduction tasks and measure recession velocities from 
emission lines using the package RVSAO (Kurtz & Mink 
1998). We flux calibrate using spectrophotometric stan- 
dards (Massey et al. 1988) and the standard Kitt Peak 
atmospheric extinction correction. For objects obtained 
in non-photometric conditions, we scale the spectra by 
the flux ratio of the observed spectroscopic and SDSS 
broadband magnitudes. We estimate that absolute flux 
calibration is accurate to ~25%. 

We measure emission line fluxes using IRAF splot 
and fitprof tasks, and find no significant offset between 
the two methods. Table 2 presents the observed line 
strengths with their measurement uncertainties. The sta- 
tistical uncertainties are formally a few percent, but the 
true error is dominated by uncertainties in the redden- 
ing correction, the stellar absorption correction, and the 
absolute flux calibration. 

For our analysis, we correct the observed emission line 
fluxes for reddening using the Balmer decrement and the 
Cardelli et al. (1989) reddening curve. We assumed an 
Ry = Av/EiB-V) = 3.1 and an intrinsic Ha/H/3 ratio 
of 2.85 (the Balmer decrement for case B recombination 
at T= lO-^K and n, - 10^ - lO'^cm--''; Osterbrock 1989). 

We find no evidence for stellar absorption in the ex- 
pected sense. While one would expect some underly- 
ing stellar absorption, we do not see Stark-broadened 
absorption in the wings of the Balmer emission lines. 
The de-reddened Balmer emission line ratios exhibit a 
7% scatter around the Osterbrock (1989) values for case 
B recombination at T= 10^ K. If we apply a constant 2 
A equivalent width correction, appropriate for the young 
4-5 Myr stellar age of the galaxies (see §3.4), the scatter 
of the Balmer line ratios around the Osterbrock (1989) 
values remains unchanged at 7%. We conclude that stel- 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which arc operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 3.— MMT spectra of our 4 new XMPGs. [O III] A4363 is 
well-detected in all four objects (see insets). The extremely weak 
[Nil] and [Sii] lines visibly indicate the extremely low abundance of 
these objects. The thumbnail images, courtesy of SDSS, are 48" on 
a side and show that all 4 XMPGs appear to be compact dwarfs. 



lar absorption is smaller than the uncertainty in the line 
ratio measurements. 

We calculate electron densities with the S II A6717 / 
S II A6731 line ratio, when present, in conjunction with a 
5 level model atom using the Mappings photoionization 
code (Sutherland fc Dopita 1993). We derive the gas- 
phase oxygen abundance following the procedure out- 
lined in Izotov et al. (2006b) within the framework of 
the classical two-zone HII- region model (Stasinska 1980). 
This procedure utilizes the electron-temperature Te cal- 
ibrations of Aller (1984) and the atomic data compiled 
by Stasinska (2005). 

The gas-phase oxygen abundance depends on line ra- 
tios and thus is independent of the uncertainties in our 
absolute fiux calibration. We propagate the errors from 
the line ratio measurement, the extinction correction, 
and the stellar absorption correction, and find that the 
relative errors in metallicities derived using the same 
method are formally < 0.07 dex. However, the absolute 
error is at least ~ 0.1 dex. Thus systematics dominate 
the errors; our metallicities are accurate at the ±0.1 dex 
level. 

2.3. Survey Efficiency 

Our initial survey of 24 candidates contains 10 metal 
poor galaxies with log(0/H) -|- 12< 8 of which 4 are 
new XMPGs. Thus our strategy is --20% efficient 
for selecting XMPGs. Spectra and thumbnail im- 
ages of the 4 new XMPGs are displayed in Figure 3 
(see also Table 3). Three of the new XMPGs (SDSS 
J142250. 72-^514516.5, SDSS J144158.32-h291434.2 and 
SDSS J225900.86-M41343.5) have lower abundance than 
SDSS J0809-I-1729. The remaining objects in our survey 
are either A stars in the Milky Way, galaxies with mod- 
est emission lines, and a few odd objects (one E-|-A, one 
possible BL Lac, one possible quasar) listed in Table 1. 

3. GALAXY PROPERTIES 

We now open our discussion to include all the 

metal poor galaxies in our survey. Including SDSS 
J0809-hl729, our survey contains 5 XMPGs and 6 metal 
poor galaxies. Table 3 summarizes the spectroscopic 
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Fig. 4. — Color-magnitude distribution of nearby BCGs (Gil 
de Paz et al. 2003) and our 11 metal poor galaxies, for wliieh we 
estimate [B — R) from SDSS photometry (Fukugita et al. 1996). 



measurements for the 11 galaxies. Because the galax- 
ies are compact objects, our spectra provide reasonable 
global estimates of their properties. We compare the 
properties of our galaxies with similar samples of blue 
compact galaxies (BCGs) and with nearby GRB hosts. 
Throughout this paper we adopt a flat, A-dominated cos- 
mology with Ho = 70 and flm = 0.3. 
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Fig. 5. — Redshift and luminosity distribution of BCGs (Kong & 
Cheng 2002; Gil de Paz et al. 2003), metal poor galaxies from 2dF 
(Papaderos et al. 2006) and SDSS (Kniazev et al. 2003), our 11 
metal poor galaxies, and the 3 nearest GRB host galaxies (Stanek 
et al. 2006). Mb is calculated for the 2dF galaxies assuming (B — 
V) = 0.5 (Papaderos et al. 2006). Mb is estimated for the SDSS 
galaxies from SDSS photometry (Kniaaev et al. 2003). 



ies at 2 ~ 0.04. Our sample of metal poor galaxies, by 
construction, spans the range 0.02 < z < 0.08. Search- 
ing a large volume of space enables the discovery of rare 
objects like XMPGs. 



3.1. Color and Redshift Distribution 

Our 11 metal poor galaxies are systematically bluer 
and more luminous than other known BCGs. Figure 4 
plots the Gil de Paz et al. (2003) sample of nearby BCGs, 
for which {B — R)q and Mb are all available. We esti- 
mate {B — R) for our galaxies from SDSS photometry 
(Fukugita et al. 1996), and shift the observed magni- 
tudes to the rest-frame using fc-corrections we calculate 
for B and R passbands as described in §2. It is clear that 
our metal poor galaxies are more than 0.5 mag bluer in 
{B - R)o than most nearby BCGs. I Zw 18 and UCM 
1612+1308 (Rego et al. 1998) are the two BCGs with 
colors comparable to our metal poor galaxies. However, 
our metal poor galaxies are systematically more luminous 
than the BCGs with similar colors. 

Our 11 metal poor galaxies are also at greater redshift 
than most known BCGs. Figure 5 plots the redshift dis- 
tribution of BCGs (Kong & Cheng 2002; Gil de Paz et al. 
2003), metal poor galaxies from 2dF (Papaderos et al. 
2006) and SDSS (Kniazev et al. 2003), our 11 metal poor 
galaxies, and the 3 nearest GRB host galaxies (Stanek 
et al. 2006). Note that we calculate Mb for the 2dF 
galaxies assuming the average (B — V) = 0.5 for that 
sample (Papaderos et al. 2006). We estimate Mb for 
the SDSS galaxies from (/ and r' photometry in Kniazev 
et al. (2003). Figure 5 shows that the vast majority of 
known BCGs have redshifts z < 0.02. The 2dF and SDSS 
surveys access fainter magnitudes than earlier BCG sam- 
ples and thus contain a number of very low luminosity 
galaxies, plus a few higher luminosity metal poor galax- 



3.2. Luminosity-Metallicity Relation 

Luminosity-metallicity relations are well determined 

for galaxies ranging from large star-forming galaxies (e.g. 
Tremonti et al. 2004) to dwarf irregulars (e.g. Richer & 
McCall 1995). The physical basis for the luminosity- 
metallicity relation is a mass-metallicity relation: low 
mass galaxies are thought to sustain less star formation 
and retain fewer metals than high mass galaxies. We now 
compare the luminosity and metallicity of our galaxies 
with samples of comparable metal-poor dwarf galaxies. 
Because there is some disagreement about the calibra- 
tion of different metallicity-estimate methods, we only 
consider galaxy samples with metallicities derived with 
the Te method. 

Figure 6 plots samples of BCGs (Kong & Cheng 2002; 
Shi et al. 2005), metal poor galaxies from 2dF (Pa- 
paderos et al. 2006) and SDSS (Kniazev et al. 2003), our 
metal poor galaxies, and four nearby GRB host galaxies 
(Stanek et al. 2006). We use metallicities for the GRB 
hosts calculated using the Te method as described in 
Kewley ct al. (2007). Figure 6 shows that our sample 
of galaxies has either 1) lower metallicity by ~0.5 dex 
than the Richer & McCall (1995) luminosity-metallicity 
relation for normal dirrs, or 2) higher luminosity by 3 - 
5 mag in Mb- 

Other samples of nearby, metal poor galaxies exhibit 
a large scatter around the luminosity-metallicity relation 
(Kunth & Ostlin 2000). Yet our galaxies appear un- 
usual because of their large systematic offset from the 
luminosity-metallicity relation. We will use population 



MMT Extremely Low Metallicity Galaxy Survey I. 



• 


Our galaxies 


I I I 1 I I I 1 I I 


□ 


GRB hosts 






I Zw 18 


- 


- o 


SBS 0335-05gW 


+ 

+ + + 


* 


SDSS (Kniazev et al) 






2dF (Papaderos et al) 


+ - ^ 4-1- + 


+ 


BCGs (Shi/Kong) 






Richer & McCall_|_ 


y ++ - 


- 


+ 


^ + + 


- 




+ •□++ + 






^* • n 






□ 










* 








• * 






• • 




% 


< < < 1 < < < 1 < < 



8.5 - 



+ 



* Our galaxies 
□ GRB hosts 

* I Zw 18 

O SBS 0335-05gW 

* SDSS (Kniazev et al) 
A 2dF (Papaderos et al) 
+ BCGs (Shi/Kong) + 

- - Richer & McCall+O.a 



V + 



7.5 - 



O 



8. 

X 

\ 

o 

o 

+ 
cv 

7. 



-10 -12 -14 -16 -18 -20 -22 

Mg [mag] 

Fig. 6. — Luminosity-mctallicity plot for BCGs (Kong & Cheng 
2002; Shi ct al. 2005), metal poor galaxies from 2dF (Papaderos 
et al. 2006) and SDSS (Kniazev et al. 2003), our metal poor galax- 
ies, and four nearby GRB host galaxies (Stanck et al. 2006). Our 
sample of metal poor galaxies fills in the region defined by the GRB 
hosts; all these galaxies have lower metallicity by ~0.5 dex than 
the Richer & McCall (1995) luminosity-metallicity relation {dashed 
line) for normal dirrs. 
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Fk;. 7. — Same as Figure 6, except we plot strong line metal- 
licitics calculated using R23 McGaugli (1991). The strong line 
metallicities are systematically +0.2 dex larger than the Te metal- 
licities, but the overall distribution remains the same: our sample 
of metal poor galaxies and the GRB hosts are systematically oflFset 
from the Richer & McCall (1995) luminosity-metallicity relation 
{dashed line), which we shift by +0.2 dex for consistency. 



synthesis models to address the evolutionary paths of 
our metal poor galaxies in a future paper. 

Remarkably, our sample of metal poor galaxies fills 
the region of the luminosity-metallicity diagram outlined 
by nearby GRB host galaxies. The only other galaxies 
with similar properties are the five metal poor 2dF and 
SDSS galaxies at z - 0.04: 2dF 169299, 2dF 115901, 
SDSSJ051902.64-h000730.0, SDSSJ104457.84+035313.2, 
and SDSSJ084030.00-F470710.2. The nearby BCG II Zw 
70 also appears to fall in the region defined by the GRB 
hosts. 

The relative distribution of galaxies in the luminosity- 
metallicity plot does not change with different metallicity 
estimators. To illustrate this point, we calculate strong 
line metallicities using R23 (McGaugh 1991) for the en- 
tire set of galaxies. The results are plotted in Figure 7. 
The R23 metallicities have an average offset of +0.2 dex 
from the Te metallicities, and so we shift the Richer & 
McCall (1995) luminosity-metallicity relation in Figure 
7 by a constant +0.2 dex for consistency. Yet the rela- 
tive distribution of galaxies remains the same: our metal 
poor galaxies (and the GRB host galaxies) maintain a 
large, systematic offset from the luminosity-metallicity 
relation defined by normal dIrrs and BGGs. 

3.3. Extinction 

Our sample of 11 metal poor galaxies suffer from very 
little internal extinction. Table 3 lists the Balmer Ha/H/? 
ratios, which average 2.8±0.2. The intrinsic Ha/H/? ratio 
is 2.86 for case B recombination at T= 10'' K and rie ~ 
10^ - 10^ cm-3 (Osterbrock 1989). Low extinction is 
expected in low metallicity galaxies, and is also observed 
in nearby GRB hosts (e.g. Kewley et al. 2007). 



3.4. Star Formation Age and Rates 

XMPGs remain a puzzle because they may be pris- 
tine galaxies undergoing their first burst of star forma- 
tion or they may contain older stellar population from 
previous episodes of star formation. Hubble Space Tele- 
scope images resolve old stellar populations in the nearest 
metal poor galaxies. I Zw 18, for example, has stars with 
ages ranging from ~ 500 Myr (Izotov & Thuan 2004) to 
^ 1 Gyr (Aloisi et al. 1999). Although we cannot esti- 
mate the age of old stellar populations (if any) in our 
metal poor galaxies, stellar population synthesis models 
provide an estimate of the age of the young stellar pop- 
ulation. 

We estimate the stellar age of otir metal poor galaxies 
from the H/3 eqtiivalent width following Schaerer & Vaeca 
(1998). Under the assumption of a Salpeter initial mass 
function and an instantaneous burst of star formation, 
the stellar ages of our metal poor galaxies are in the 
range 4-5 Myr (see Table 3). 

Star formation rates (SFRs) are more difficult to esti- 
mate because traditional calibrations are not applicable 
to extremely metal poor objects. We use the calibration 
derived by Kewley et al. (2007) based on the stellar popu- 
lation synthesis models of Bicker & Ftitze-v. Alvensleben 
(2005). Our metal poor galaxies then have SFRs in the 
range 0.1 - 0.5 M0 yr^^ (Table 3). Both the stellar ages 
and SFRs of our metal poor galaxies are remarkably sim- 
ilar to those found for nearby GRB hosts (Kewley et al. 
2007). 

4. SUPERNOVA RATES AND THE GRB CONNECTION 

We next ask how frequently GRBs might occur in our 

sample of metal poor galaxies, if GRBs are indeed associ- 
ated with core-collapse supernova in metal poor galaxies. 
We start by estimating the number of massive O and B 
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stars in our sample of 11 galaxies, and then estimate the 
rate of core-collapse supernovac. 

Under the assumption of an ionization-bounded neb- 
ula, the H/3 line luminosity provides an estimate of the 
ionizing flux present in the galaxies. The metal poor 
galaxies suffer from very little extinction, as measured 
by their Ha/H/? ratios; thus we make no correction for 
internal extinction. We follow Schaerer & Vacca (1998) 
and convert the observed H/3 ionizing flux to an equiva- 
lent number of O stars. We select the "equivalent 07V 
to O star" ratio r]o based on the stellar age and metal- 
licity of the galaxy; 770 values range from 0.25 - 0.5. The 
minimum O star mass is ^^13. 3 Mq at our metallicities 
(Vacca 1994). However, core-collapse supernovae will 
result from B stars with masses as low as ~8 Mq. A 
Salpeter initial mass function has the same number of O 
stars with 13 - 120 M© as B stars with 8 - 13 Mq. Thus, 
the total number of core-collapse supernova progenitors 
is roughly twice the Schaerer & Vacca (1998) number of 
equivalent O stars. We list our estimate of the total num- 
ber of core-collapse supernova progenitors in the column 
n(OB) (Table 3). On average, there arc ^^25,000 poten- 
tial core-collapse supernova progenitors per galaxy in our 
sample. This number may appear relatively small, but 
the metal poor galaxies arc dwarfs, not massive galaxies. 

We estimate core-coUapsc supernova rates in our metal 
poor galaxies by assuming supernovae occur uniformly 
over the lifetime of the longest-lived progenitor. The 
Schaller et al. (1992) stellar evolution track for a 9 Mq 
star with Z = 0.001 has a hfetime of 30 Myr. Thus 
the average rate of core-collapse supernovae in our metal 
poor galaxies is ~10~^ yr""'^, with an uncertainty of a 
factor of a few. In other words, a couple thousand such 
metal poor galaxies must be monitored to witness one 
core-collapse supernova per year. A typical spiral galaxy, 
in comparison, has a typical SFR of ~10 M© yr~-^ (e.g. 
Brinchmann et al. 2004) and a core-coUapsc supernova 
rate ~100 times larger than the metal poor galaxies, as- 
suming a constant SFR and Salpeter initial mass func- 
tion. 

The ratio of GRBs to core-collapse supernovae events 
depends on how narrowly beamed GRBs are. If 1 out 
of 100 core-collapsc supernovae appear as GRBs, then a 
sample of a ^^10"'' metal poor galaxies must be monitored 
to witness one GRB per year. There is currently little 
constraint on the space density of metal poor galaxies, in 
part because they are such low-luminosity systems. Sur- 
veys reaching faint magnitude limits may uncover large 
numbers of new metal poor galaxies. However, it may 
be difhcult to detect a supernova coincident on the high 
surface-brightness core of these galaxies. Figure 3 shows 
that the four XMPGs, for example, are compact sys- 
tems. Although metal poor galaxies are potential GRB 
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hosts, a survey to find GRBs by monitoring a sample 
metal poor galaxies appears impractical currently. Fu- 
ture deep imaging surveys, such as Pan-STARRS and the 
Large Synoptic Survey Telescope (LSST), may be able to 
detect supernovae in these low-luminosity galaxies. 

5. CONCLUSIONS 

We have designed a successful strategy to find new 
metal poor galaxies. We calculate the expected col- 
ors of metal poor galaxies at different redshifts based 

on the spectrum of the newly discovered XMPG, SDSS 
J0809-hl729 (Kewley et al. 2007). We observed an ini- 
tial sample of 24 candidates with the MMT telescope and 
find 4 new XMPGs with log(0/H) + 12< 7.65, a ~20% 
selection efficiency. 

Our full set of 11 metal poor galaxies are systemat- 
ically bluer and more luminous than comparable sam- 
ples of BCGs. Our galaxies are also systematically more 
metal poor by 0.5 dex (or more luminous by 3 - 5 mag) 
than samples of BCGs, dirrs, and other metal-poor dwarf 
galaxies. Remarkably, our galaxies share the same region 
of the luminosity-metallicity diagram with nearby GRB 
hosts. The similarity in extinction, stellar age, and star 
formation rates suggests that our metal poor galaxies are 
potential hosts for GRBs. 

We estimate an average core-collapse supernova rate 
~10~^ yr~^ in our metal poor galaxies. This estimate 
comes from an estimate of the number of O and B stars in 
the galaxies. If GRBs arc indeed linked to core-collapse 
supernova in metal poor galaxies, future surveys such as 
Pan-STARRS or LSST may be able to find GRBs by 
monitoring a large sample of metal poor galaxies. 

The success of our XMPG selection strategy allows 
us to expand our survey. For example, using the new 
XMPGs (Figure 3) as additional templates, we identify 
a total of 335 XMPG candidates in SDSS Data Release 
5. Spectroscopic observations of these XMPG candidates 
are underway. 
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TABLE 1 
XMPG CANDIDATES 





ID 


9 

(mag) 


cz 

/I — 1 \ 
(km s ^) 


Comment 


SDSS 




18.70 






SDSS 




20.17 


23340 




SDSS 


T12SQ44 ^8+14^1(112 8 


19.77 


21534 




SDSS 


T1246S8 82+'^5mi5 1 


20.36 


19522 


lilt- LClfl v^VJl. &C*ilCUrt. V 


SDSS 


J124709. 24+325118. 6 


19.10 


27450 


era 1 a"5rv 


SDSS 


T1 3S424 53+5Q2n57 fl 


20.55 


21890 


lll^Ldl VJKJXjl. gdilCUA. V 


SDSS 


T1 35fi41 fi4+fi54748 8 


19.72 


10742 




SDSS 


T14n4'HQ 28+542136 Q 


19.98 


318 




SDSS 


T1 41 333 54+463234 D 


19.43 


90 




SDSS 


1142250 72+51 451 fi 5 


20.22 


11654 




SDSS 


T1 43345 QQ-n25fin2 2 


19.69 


65 




SDSS 


T1 441 58 32+2Q1 434 2 


20.13 


13741 




SDSS 


J145621.69+503523.0 


20.49 




BL Lac? 


SDSS 


J150316.52+111056.9 


19.41 


23405 


metal poor galaxy 


SDSS 


J150535.89+314639.4 


20.55 


15817 


galaxy 


SDSS 


J151221.08+054911.2 


20.15 


24025 


metal poor galaxy 


SDSS 


J152802.62+240425.6 


18.65 


-6 


A-star 


SDSS 


J154742.23-005554.2 


20.57 


16594 


E+A galaxy 


SDSS 


J160238.71+444923.8 


19.65 


12350 


galaxy 


SDSS 


J165835. 08+192415.3 


18.93 


-23 


A-star 


SDSS 


,1172955.61+534338.8 


19.55 


24308 


metal poor galaxy 


SDSS 


J211613. 97-000851. 3 


20.39 


11916 


galaxy 


SDSS 


J221912. 56+140602.8 


20.41 


36 


A-star 


SDSS 


J225900.86+141343.5 


19.09 


8918 


metal poor galaxy 
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TABLE 2 
OBSERVED LINE INTENSITIES 



F(Ao Ion)=' 


J120955.67 


J123944.58 


J124638.82 


J133424.53 


J142250.72 


J144158.32 


J150316.52 


J151221.08 


J172955.61 


J22E 


+142155.9 


+145612.8 


+350115.1 


+592057.0 


+514516.5 


+291434.2 


+111056.9 


+054911.2 


+534338.8 


+14 



3727 [On] 
3798 HID 
3835 H9 
3868 [Nciii] 
3889 He I + 
3968 [Nciii] 
4101 H(5 
4340 H7 
4363 [Oiii] 
4471 He I 
4686 He 11 
4861 H/3 
4959 [Oiii] 
5007 [Oiii] 
5876 He l 
6563 Ha 
6584 [Nil] 
6678 He i 
6717 [Sii] 
6731 [Sii] 
EW(H^) 



8.20 ± 0.22 



H8 
+ H7 



1, 


.75 




0.20 


0, 


.66 


± 


0.08 


1, 


.07 


± 


0.11 


1, 


.40 


± 


0.11 


2, 


.70 


± 


0.10 


0, 


.30 


± 


0.04 



5.50 ± 0.11 
5.50 ± 0.24 

16.3 ± 0.2 
0.46 ± 0.06 

14.8 ± 0.2 



1.40 ± 0.26 
0.50 ± 0.08 
54.4 ± 3.0 



11.8 ± 0.2 
0.84 ± 0.10 
0.29 ± 0.08 
3.38 ± 0.14 
2.12 ± 0.17 
2.50 ± 0.11 
2.50 ± 0.12 
4.30 ± 0.11 
1.10 ± 0.13 



10.4 ±0.1 
15.1 ± 0.4 
45.0 ± 0.5 
2.60 ± 0.11 
30.7 ± 0.4 



1.00 ± 0.16 
1.00 ± 0.15 
91.2 ± 5.0 



16.4 ± 0.4 
0.49 ± 0.06 
0.34 ± 0.06 
4.09 ± 0.13 



1.86 : 

2.39 

2.40 

4.40 

1.00 

0.36 

0.19 : 



0.21 
0.12 
0.23 
0.10 
0.11 
0.06 
0.05 



± 
± 
± 
± 
± 
± 
± 

10.3 ± 0.1 
14.7 ± 0.1 

44.6 ± 0.5 
1.18 ± 0.11 
31.5 ± 0.3 
0.70 ± 0.08 

1.40 ± 0.11 
0.90 ± 0.13 

70.7 ± 3.9 



22.5 ± 0.3 
0.62 ± 0.07 
0.55 ± 0.07 
5.62 ± 0.11 
2.22 ± 0.11 
3.84 ± 0.14 
3.60 ± 0.18 
6.10 ± 0.10 
1.10 ± 0.10 
0.32 ± 0.04 
0.23 ± 0.04 
14.1 ± 0.1 
20.8 ± 0.2 
61.8 ± 0.7 
1.73 ± 0.12 
40.3 ± 0.7 
1.30 ± 0.21 
0.39 ± 0.06 
2.20 ± 0.11 
1.30 ± 0.12 
72.5 ± 4.0 



7.80 ± 0.14 
0.35 ± 0.07 
0.91 ± 0.10 
6.11 ± 0.27 
3.14 ± 0.11 
4.78 ± 0.28 
3.80 ± 0.11 
8.30 ± 0.11 
2.00 ± 0.20 
0.73 ± 0.09 
0.73 ± 0.09 
18.3 ± 0.2 
21.8 ± 0.2 
65.5 ± 0.7 
2.08 ± 0.32 
46.2 ± 0.5 



0.70 ± 0.08 
0.50 ± 0.09 
95.9 ± 5.3 



5.70 ± 0.13 
0.27 ± 0.05 
0.45 ± 0.07 
1.47 ± 0.12 
0.84 ± 0.10 
0.89 ± 0.09 
0.90 ± 0.11 

I. 90 ± 0.11 
0.40 ± 0.07 
0.29 ± 0.04 

4.00 ± 0.10 
4.50 ± 0.21 

13.6 ± 0.2 
0.37 ± 0.05 

II. 3 ± 0.1 
0.20 ± 0.03 

0.50 ± 0.09 
0.30 ± 0.06 

76.7 ± 4.2 



43.8 ± 0.5 
0.72 ± 0.09 
1.35 ± 0.12 

14.5 ± 0.3 

5.43 ± 0.11 
8.67 ± 0.15 
8.50 ± 0.19 

15.9 ± 0.3 
3.50 ± 0.14 
1.31 ± 0.11 

37.9 ± 0.4 
62.8 ± 0.6 
190.0 ± 1.9 
3.73 ± 0.12 
102.0 ± 1.2 
1.90 ± 0.11 
0.94 ± 0.25 
3.70 ± 0.12 
3.20 ± 0.11 
80 ± 4.4 



9.00 ± 0.12 
0.36 ± 0.10 
0.35 ± 0.06 
2.98 ± 0.23 
1.44 ± 0.11 

I. 90 ± 0.11 
2.30 ± 0.13 
4.00 ± 0.11 
0.70 ± 0.09 

0.37 ± 0.05 
9.80 ± 0.15 

II. 9 ± 0.2 
35.2 ± 0.4 
1.03 ± 0.11 
25.8 ± 0.3 
0.40 ± 0.05 

0.90 ± 0.10 

61.6 ± 3.4 



29.9 ± 0.5 
1.25 ± 0.11 
1.38 ± 0.11 

13.5 ± 0.4 
5.09 ± 0.29 
8.21 ± 0.11 
6.90 ± 0.14 

12.4 ± 0.2 
2.90 ± 0.19 

0.41 ± 0.09 
29.2 ± 0.4 

62.5 ± 0.7 
188.0 ± 2.2 
3.76 ± 0.13 
94.2 ± 1.0 
1.40 ± 0.11 
0.83 ± 0.09 
3.20 ± 0.14 
2.20 ± 0.12 

158 ± 8.7 



15 
1.2 
1.3 
7.0 
5.5 
6.8 
7.9 
15 
2.9 
0.6 
0.6 
33 
33 
10( 
2.8 
91 
0.4 
0.4 
1.4 
1.0 
13 



Tlux in units of X 10 
calibration. 



ergs s cm 



Uncertainties are statistical errors only, and do not include the errors due to reddening, stellar absorption, and absc 



TABLE 3 
OUR METAL POOR GALAXIES 





ID 


cz 


Mb 




log(0/H) + 12 


Age 


Ha/H/3 


L(H^) 


SFR (Te) 


N(OB) 


SN rate 






(km s^^) 


(mag) 


(K) 




(Myr) 




(ergs 1) 


(Mo yr-i) 




(10-" yr-l) 


SDSS 


J225900.86+141343.51 


8918 


-16.4 


18600 


7.37 


4.0 


2.75 


1.80x10""' 


0.21 


13900 


4 


SDSS 


J142250.72+514516.49 


11654 


-15.9 


19100 


7.41 


4.0 


2.52 


8.46x10'^'' 


0.09 


6500 


2 


SDSS 


J144158.32+291434.22 


13741 


-16.3 


19400 


7.47 


4.5 


2.80 


1.08x10*" 


0.13 


10400 


3 


SDSS 


J080840.85+172856.48'' 


13233 


-17.1 


19000 


7.48 


4.5 


2.64 


1.94x10'"' 


0.18 


18700 


6 


SDSS 


J123944.58+145612.80 


21534 


-17.7 


17000 


7.65 


4.5 


2.95 


3.69x10'"' 


0.54 


37400 


12 


SDSS 


J151221. 08+054911. 21 


24025 


-17.6 


15300 


7.68 


5.5 


2.64 


2.63x10'"' 


0.32 


40600 


14 


SDSS 


J120955.67+142155.91 


23340 


-17.4 


15000 


7.71 


5.5 


2.72 


2.26x10'"' 


0.29 


34900 


12 


SDSS 


J124638.82+350115.11 


19522 


-16.8 


16200 


7.75 


5.0 


3.06 


1.61x10"" 


0.26 


24900 


8 


SDSS 


J150316.52+111056.93 


23405 


-18.3 


14800 


7.85 


4.5 


2.69 


6.31x10"" 


0.80 


65800 


22 


SDSS 


J133424.53+592057.04 


21890 


-16.8 


14500 


7.86 


5.0 


2.86 


1.73x10"" 


0.23 


27800 


9 


SDSS 


J172955.61+534338.80 


24308 


-18.3 


13900 


8.00 


4.0 


3.23 


9.46x10"" 


1.98 


60800 


20 
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